Detection of selective sweeps in structured populations : a comparison of recent methods

This work was supported by the Marie-Curie Initial Training Network INTERCROSSING (European Commission FP7). OEG was further supported by the MASTS pooling initiative (The Marine Alliance for Science and Technology for Scotland). Date of Acceptance: 25/08/2015Identifying genomic regions targeted by positive selection has been a longstanding interest of evolutionary biologists. This objective was difficult to achieve until the recent emergence of Next Generation Sequencing, which is fostering the development of large-scale catalogs of genetic variation for increasing number of species. Several statistical methods have been recently developed to analyze these rich datasets but there is still a poor understanding of the conditions under which these methods produce reliable results. This study aims at filling this gap by assessing the performance of genome-scan methods that consider explicitly the physical linkage among SNPs surrounding a selected variant. Our study compares the performance of seven recent methods for the detection of selective sweeps (iHS, nSL, EHHST, xp-EHH, XP-EHHST, XPCLR and hapFLK). We use an individual-based simulation approach to investigate the power and accuracy of these methods under a wide range of population models under both hard and soft sweeps. Our results indicate that XPCLR and hapFLK perform best and can detect soft sweeps under simple population structure scenarios if migration rate is low. All methods perform poorly with moderate to high migration rates, or with weak selection and very poorly under a hierarchical population structure. Finally, no single method is able to detect both starting and nearly completed selective sweeps. However, combining several methods (XPCLR or hapFLK with iHS or nSL) can greatly increase the power to pinpoint the selected region.PostprintPeer reviewe

A complex pattern of post‐divergence expansion, contraction, introgression and asynchronous responses to Pleistocene climate changes in two Dipelta sister species from western China

The well-known vicariance and dispersal models dominate in understanding the allopatric pattern for related species and presume the simultaneous occurrence of speciation and biogeographic events. However, the formation of allopatry may postdate the species divergence. We examined this hypothesis using DNA sequence data from 3 chloroplast fragments and 5 nuclear loci of Dipelta floribunda and D. yunnanensis, two shrub species with the circum Sichuan Basin distribution, combining the climatic niche modeling approach. The best-fit model supported by the approximate Bayesian computation (ABC) analysis indicated that, D. floribunda and D. yunnanensis diverged during the mid-Pleistocene period, consistent with the largest glacial period in the Qinghai-Tibet Plateau (QTP). The historically inter-specific gene flow was identified but seemed to have ceased after the last interglacial period (LIG), when the range of D. floribunda moved northward from the south of the Sichuan Basin. Further, populations of D. floribunda had expanded obviously in the north of the Sichuan Basin after the last glacial maximum (LGM). Relatively, the range of D. yunnanensis expanded before the LGM, reduced during the post-LGM especially in the north of the Sichuan Basin, reflecting the asynchronous responses of related species to the contemporary climate changes. Our results suggested that complex topography should be considered in understanding the distributional patterns even for closely related species and their demographic responses

Fluctuating selection models and Mcdonald-Kreitman type analyses

It is likely that the strength of selection acting upon a mutation varies through time due to changes in the environment. However, most population genetic theory assumes that the strength of selection remains constant. Here we investigate the consequences of fluctuating selection pressures on the quantification of adaptive evolution using McDonald-Kreitman (MK) style approaches. In agreement with previous work, we show that fluctuating selection can generate evidence of adaptive evolution even when the expected strength of selection on a mutation is zero. However, we also find that the mutations, which contribute to both polymorphism and divergence tend, on average, to be positively selected during their lifetime, under fluctuating selection models. This is because mutations that fluctuate, by chance, to positive selected values, tend to reach higher frequencies in the population than those that fluctuate towards negative values. Hence the evidence of positive adaptive evolution detected under a fluctuating selection model by MK type approaches is genuine since fixed mutations tend to be advantageous on average during their lifetime. Never-the-less we show that methods tend to underestimate the rate of adaptive evolution when selection fluctuates

Analysis of IL2/IL21 Gene Variants in Cholestatic Liver Diseases Reveals an Association with Primary Sclerosing Cholangitis

Background/Aims: The chromosome 4q27 region harboring IL2 and IL21 is an established risk locus for ulcerative colitis (UC) and various other autoimmune diseases. Considering the strong coincidence of primary sclerosing cholangitis (PSC) with UC and the increased frequency of other autoimmune disorders in patients with primary biliary cirrhosis (PBC), we investigated whether genetic variation in the IL2/IL21 region may also modulate the susceptibility to these two rare cholestatic liver diseases. Methods: Four strongly UC-associated single nucleotide polymorphisms (SNPs) within the KIAA1109/TENR/IL2/IL21 linkage disequilibrium block were genotyped in 124 PBC and 41 PSC patients. Control allele frequencies from 1,487 healthy, unrelated Caucasians were available from a previous UC association study. Results: The minor alleles of all four markers were associated with a decreased susceptibility to PSC (rs13151961: p = 0.013, odds ratio (OR) 0.34; rs13119723: p = 0.023, OR 0.40; rs6822844: p = 0.031, OR 0.41; rs6840978: p = 0.043, OR 0.46). Moreover, a haplotype consisting of the four minor alleles also had a protective effect on PSC susceptibility (p = 0.0084, OR 0.28). A haplotype of the four major alleles was independently associated with PSC when excluding the patients with concomitant inflammatory bowel disease (p = 0.033, OR 4.18). Conclusion: The IL2/IL21 region may be one of the highly suggestive but so far rarely identified shared susceptibility loci for PSC and UC. Copyright (C) 2011 S. Karger AG, Base

Genome-Wide Analysis of Natural Selection on Human Cis-Elements

Background: It has been speculated that the polymorphisms in the non-coding portion of the human genome underlie much of the phenotypic variability among humans and between humans and other primates. If so, these genomic regions may be undergoing rapid evolutionary change, due in part to natural selection. However, the non-coding region is a heterogeneous mix of functional and non-functional regions. Furthermore, the functional regions are comprised of a variety of different types of elements, each under potentially different selection regimes. Findings and Conclusions: Using the HapMap and Perlegen polymorphism data that map to a stringent set of putative binding sites in human proximal promoters, we apply the Derived Allele Frequency distribution test of neutrality to provide evidence that many human-specific and primate-specific binding sites are likely evolving under positive selection. We also discuss inherent limitations of publicly available human SNP datasets that complicate the inference of selection pressures. Finally, we show that the genes whose proximal binding sites contain high frequency derived alleles are enriched for positive regulation of protein metabolism and developmental processes. Thus our genome-scale investigation provide

Direct Integration and Non-Perturbative Effects in Matrix Models

We show how direct integration can be used to solve the closed amplitudes of multi-cut matrix models with polynomial potentials. In the case of the cubic matrix model, we give explicit expressions for the ring of non-holomorphic modular objects that are needed to express all closed matrix model amplitudes. This allows us to integrate the holomorphic anomaly equation up to holomorphic modular terms that we fix by the gap condition up to genus four. There is an one-dimensional submanifold of the moduli space in which the spectral curve becomes the Seiberg--Witten curve and the ring reduces to the non-holomorphic modular ring of the group $\Gamma(2)$. On that submanifold, the gap conditions completely fix the holomorphic ambiguity and the model can be solved explicitly to very high genus. We use these results to make precision tests of the connection between the large order behavior of the 1/N expansion and non-perturbative effects due to instantons. Finally, we argue that a full understanding of the large genus asymptotics in the multi-cut case requires a new class of non-perturbative sectors in the matrix model.Comment: 51 pages, 8 figure

Strong signature of natural selection within an FHIT intron implicated in prostate cancer risk

Previously, a candidate gene linkage approach on brother pairs affected with prostate cancer identified a locus of prostate cancer susceptibility at D3S1234 within the fragile histidine triad gene (FHIT), a tumor suppressor that induces apoptosis. Subsequent association tests on 16 SNPs spanning approximately 381 kb surrounding D3S1234 in Americans of European descent revealed significant evidence of association for a single SNP within intron 5 of FHIT. In the current study, resequencing and genotyping within a 28.5 kb region surrounding this SNP further delineated the association with prostate cancer risk to a 15 kb region. Multiple SNPs in sequences under evolutionary constraint within intron 5 of FHIT defined several related haplotypes with an increased risk of prostate cancer in European-Americans. Strong associations were detected for a risk haplotype defined by SNPs 138543, 142413, and 152494 in all cases (Pearson's χ2 = 12.34, df 1, P = 0.00045) and for the homozygous risk haplotype defined by SNPs 144716, 142413, and 148444 in cases that shared 2 alleles identical by descent with their affected brothers (Pearson's χ2 = 11.50, df 1, P = 0.00070). In addition to highly conserved sequences encompassing SNPs 148444 and 152413, population studies revealed strong signatures of natural selection for a 1 kb window covering the SNP 144716 in two human populations, the European American (π = 0.0072, Tajima's D= 3.31, 14 SNPs) and the Japanese (π = 0.0049, Fay & Wu's H = 8.05, 14 SNPs), as well as in chimpanzees (Fay & Wu's H = 8.62, 12 SNPs). These results strongly support the involvement of the FHIT intronic region in an increased risk of prostate cancer. © 2008 Ding et al

Incipient Balancing Selection through Adaptive Loss of Aquaporins in Natural Saccharomyces cerevisiae Populations

A major goal in evolutionary biology is to understand how adaptive evolution has influenced natural variation, but identifying loci subject to positive selection has been a challenge. Here we present the adaptive loss of a pair of paralogous genes in specific Saccharomyces cerevisiae subpopulations. We mapped natural variation in freeze-thaw tolerance to two water transporters, AQY1 and AQY2, previously implicated in freeze-thaw survival. However, whereas freeze-thaw–tolerant strains harbor functional aquaporin genes, the set of sensitive strains lost aquaporin function at least 6 independent times. Several genomic signatures at AQY1 and/or AQY2 reveal low variation surrounding these loci within strains of the same haplotype, but high variation between strain groups. This is consistent with recent adaptive loss of aquaporins in subgroups of strains, leading to incipient balancing selection. We show that, although aquaporins are critical for surviving freeze-thaw stress, loss of both genes provides a major fitness advantage on high-sugar substrates common to many strains' natural niche. Strikingly, strains with non-functional alleles have also lost the ancestral requirement for aquaporins during spore formation. Thus, the antagonistic effect of aquaporin function—providing an advantage in freeze-thaw tolerance but a fitness defect for growth in high-sugar environments—contributes to the maintenance of both functional and nonfunctional alleles in S. cerevisiae. This work also shows that gene loss through multiple missense and nonsense mutations, hallmarks of pseudogenization presumed to emerge after loss of constraint, can arise through positive selection